Perturbative criteria for Anderson localization in long-ranged ID tight-binding models 
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We develop an alternative scaling approach to determine the criteria for Anderson localization in 
one-dimensional tight-binding models with random site energies having a bandwidth that decays as 
a power law in space, Hij oc \i — At the first order in perturbation theory the scale dependence 

of the exchange-narrowed energy of the disorder is compared to the energy level spacing of the ideal 
system to establish whether or not the disorder has a perturbative effect on the Bloch states. We find 
that at a = 1, the perturbative condition is satisfied and for sufficiently weak disorder strength all 
states are extended. For a > 1, all states are localized for arbitrary disorder strength, in agreement 
with the earlier renormalization group treatment by Levitov. 

PACS numbers: 63.22.-(-m,63.50.-fx,71.23.-k,72.15.Rn,73.20.Mf 



I. INTRODUCTION 

The subject of strong localization of electronic states 
was pioneered by early efforts of P. W. Anderson (1958)^, 
revealing the existence of a phase transition from a metal- 
lic phase to an insulating phase and established a re- 
markable set of criteria or bounds on the transition in 
terms of the relative disorder strength. In one-dimension 
(ID), it can be rigorously shown that there is a strict ab- 
sence of extended states for an arbitrarily finite amount 
of disorder. This was first stated by Mott and Twose 
(1961) in a heuristic manner^ and later made more rig- 
orous by others^. However, the validity of this result 
only encompasses models with uncorrelated disorder and 
nearest- neighbor hopping terms. 

Recent efforts have studied particular classes of Hamil- 
tonians containing long-ranged interactions that exhibit 
delocalization criticality in lower dimensionsiii^i^ii^. The 
power-law random-banded matrix (PRBM) ensemble, in- 
troduced by Mirlin et al,* is a large N x N random ma- 
trix whose entries Hij decrease in a power-law fashion 
Hij (X aij \i — j\ with a random hopping coupling con- 
stant aij . Mapping the problem onto a nonlinear cr-model 
with a non-local interaction, they found a transition from 
localized to extended states at a = 1. Levitov, using a 
renormalization group approach studied a similar system 
with disordered site energies and non-random hopping, 
also found that for a > 1, all states are localized^, with 
a = 1 as the critical value. 

In this article, we adopt a perturbative approach that 
is similar in spirit to the original arguments of Anderson, 
in that a breakdown of perturbation theory in the weak 
disorder limit indicates that the ground state of the dis- 
ordered system is very different from the ideal crystal and 
the two ground-states cannot be approached by adiabat- 
ically turning the disorder strength on or off. We modify 
a scaling argument originally discussed in the reference^ 
for short-ranged systems. We consider the weak disor- 
der regime and the convergence criteria that determine 
whether the disorder is allowed to be a weak perturba- 
tive effect in the thermodynamic limit is approached but 
not taken explicitly. If this criterion is rendered false as 



N increases then all of the eigenstates are localized. We 
also provide an estimate of the scaling behavior of lo- 
calization length as function of the disorder strength in 
the regime where the disorder is moderately strong and 
nonperturbative. 



II. THE NON-DISORDERED DISPERSION 

Consider the tight-binding model with one orbital per 
site, for N sites with lattice spacing a, that is governed 
by the following Schrodinger equation. 



(2.1) 



where ipi are the electron amplitudes at integer site values 
I and the summation is taken over all unique pairs of 
lattice sites z. The effective width of the spectrum has 
the power law spatial dependence. 
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where J is a constant that depends on the details of the 
atomic orbital overlap matrix elements and the hopping 
exponent is a positive integer a = 1, 2, 3, ... . It follows 
that for translationally invariant systems with periodic 
boundary conditions we can make use of the Bloch theo- 
rem, yielding a solution to Eq. (|2.ip of the form 
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exp(iK • Ri) 



(2.3) 



where we have introduced the reciprocal lattice vector 
K = 2TTk/Na, in terms of the wave vector k and Ri = la. 
Next, by substituting Eqs. f^ and into Eq. (PT|) . 
we have the following expression for the band dispersion: 
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for a mean site energy value, {eq). Eq. ()2.4p can be eval- 
uated explicitly in terms of polylogarithms as shown in 
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the referenceiii. Using the definition: 

oc ^ 
k=l 



(2.5) 



we arrive at a more compact expression, 

^^K,. = (eo> - J{Li^[e'^] + Lz„[e-*^]) (2.6) 

It important to note that the fundamental hmitation of 
this approach is that one cannot apply the identity (|2.5p 
in the case fractional values of a. Therefore, for the en- 
tirety of this paper we only consider integer values of a 
and we concede that a critical value alpha might be a non- 
integer value however, we can estimate ranges rounded 
to the nearest integer. 

For specific cases of a one can make use of the following 
exactly summable series^i, 



a = 1,5: 
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~ 90 ^ 12 T2" ~ 48 
This results in the following dispersion relations: 
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III. PERTURBATIVE CONVERGENCE IN THE 
PRESENCE OF DISORDER 

A. The model 



where we have defined Di = £/ — {eq) as the random 
deviation of the site energies from the mean (eo)- Di is 
a random variable, distributed according to a standard 
white-noise, Gaussian distribution: 



V{Di) = (27rcr2)-i/2gxp(_i:)2/2a2) 



(3.2) 



that is centered at zero {Di) 
tion a. 



0. with a standard devia- 



B. Scaling the perturbative limit 

The scaling procedure requires two major compo- 
nents. First we must determine the energy level spac- 
ing of unperturbed, ideal system, which is dominated by 
the long wavelength behavior of the dispersion relation, 
lim ^ , that we term [N) . This quantity solely 

depends on the mathematical form of the dispersion re- 
lations, which ultimately depend on the range of the in- 
teractions. Second, we must compute the statistical fluc- 
tuations of the disorder-induced energy shifts at the first 
order in perturbation theory, while assuming that any 
divergence will occur at the lowest order. 

We begin by determining the fluctuations of the first 
order energy shifts. We will examine the diagonal and 
off-diagonal contributions separately. Since the site ener- 
gies still reside on an underlying lattice we can transform 
Eq. ()3.1|) to the K-representation by multiplying both 
sides by (j2.3p and summing over the site numbers, /: 



{E - El,M{K) = 5: VkK'V'(K') (3.3) 

K' 

where we have defined the following, 

V'(K) = -^5:^,e^K/a (3^4) 



Having established the dispersion relations for the ideal 
system, let us now turn our attention the disordered sys- 
tem, where the site energies take random values from 
an uncorrelated Gaussian distribution. We shall follow 
precisely, the formalism and notation of F. Dominguez- 
Adame and V.A. Malyshev in the reference^. The equa- 
tions of motion result from the following modification to 
Eq.lEII), 



Notice that Vkk' are simply the matrix elements, of the 
random potential, or site energies in the unperturbed ba- 
sis. The diagonal elements Fkk represent a first order 
perturbative energy shift of the 'K.th state, 

N-l 
/=0 



Dii^i + (£0) - 51 = EiPi (3.1) 
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which solely depends on the particular quenched realiza- 
tion of the disorder. The square-root fluctuations are 
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FIG. 1: A comparison of the scale dependence of the energy 
level spacings with the exchange-narrowed fluctuations of the 
disordered site-energies. Evidently, at a = 1 the Bloch states 
of the ideal system satisfy the perturbative condition, increas- 
ing faster with increasing A'' than the critical separatrix. 



given by, 
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Evidently, the diagonal contributions result in an inho- 
mogeneous broadening of the quasi-particle energy level, 
for given ensemble of quenched disorder. 

On the other hand, the off-diagonal components are 
generally associated with overlapping, or inter-site scat- 
tering of the quasi-particle states. Thus, playing an im- 
portant role in localization. By utilizing Eq. (|3.5p . one 
can see that 
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In addition, this represents the typical fluctuation, for a 
given lattice size. Notice that these fluctuations vanish 
in the limit of a large system size N ^ oo. This effect is 
known as "self-averaging" or "exchange-narrowing" . 

Most importantly, our test for localization is the per- 
turbative condition below: 
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(3.9) 



The left hand side of side of Ea. (|3.9p has the following 



scale dependence for a given value of a, 

5E^ ^ ~ -2Jln(K) cx Jln(7V) (3.10) 
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(3.12) 



Let us now examine the energy level spacings and if they 
satisfy the condition p.9p . for increasing system size N, 



lim SE<^ , » 
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(3.13) 
(3.14) 
(3.15) 



Apparently, for a — 1 the perturbative condition is satis- 
fied with increasing N, therefore the mixing of the Bloch 
states is weak enough such that the quasi-particle states 
are expected to remain extended over the entire system 
length. However, for a = 2, 4 the disorder has a non- 
perturbative effect and the perturbative condition (|3.9|) 
is violated such that the Bloch states are strongly mixed 
and the quantity decreases slower than SE^ ^ ^^^d 

5E^ 4 signifying that the disordered eigenfunctions can- 
not be perturbed from the extended Bloch states and 
thus, all of the eigenstates are localized. A useful illus- 
tration is shown in Fig[l] where the scale dependence of 
the energy level spacings for each a case is plotted on a 
logarithmic scale. Clearly, as the range of the interac- 
tions is decreased with increasing a, the eigenstates are 
more strongly localized. 



IV. SCALING AT THE ONSET OF 
LOCALIZATION 

A. The localization length 

The threshold at which the condition p.9p fails indi- 
cates the onset of localization, where localization length 
is smaller than the size of the lattice. This clearly occurs 
in the non-perturbative regime, where strength of the dis- 
order is strong but finite and intermediate between the 
weak and strong disorder regimes. Therefore, one can 
estimate how the localization length scales as a function 
of the ratio J/cr, by solving for N = m the equation 
below. 
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The solution to Eq. (l4.ip . can be understood graphically 
in FiglT] as the intersection of each of the a — 1,2,4 
curves with the separatrix. 
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For the a = 1 case, the locahzation length solves the 
equation below, 



ZTT 



(4.2) 



It follows that for sufficiently weak disorder, the ratio 
ct/J ^ I, therefore Eq. (j4.2p reduces to the cubic equa- 
tion. 
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(4.3) 



which, admits nonreal, complex solutions and as ex- 
pected the localization length should not scale with a/ J 
given that all the states are extended. On the other hand, 
in the opposite limit, a/J^ 1, the solution involves an 
essential singularity, also admitting complex solutions, 
therefore the proper scaling of the localization occurs in 
the highly nonperturbative regime, where our Eq. ()4.2|) 
must be determined numerically. This breakdown of the 
localization length scaling in both asymptotic limits sug- 
gests there is a phase transition connecting the strong 
and weak disorder regimes. 

Next for the strictly localized cases a = 2, 4, for shorter 
ranged hopping, the scaling of the localization length be- 
comes, 
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(4.4) 
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disorder, for instance, with power law decay. Instead 
of a white noise distribution of the random site ener- 
gies, consider a scenario where the power spectrum of 
the noise has S[k) cx l/Zc^, recently discussed as possibly 
occurring in conducting DNA wires and in conducting ID 
systems^^. Evidently, the autocorrelation function of the 
matrix elements of Vkk' as designated by Eq. (j3.5p would 
not scale as l/vW, rather the fluctuations of the disorder 
would decrease faster, depending on the value of j3. This 
implies that for longer ranged correlated site energies, 
the eigenstates can more easily satisfy the perturbative 
condition given by Eq. (|3.9p . and thus longer ranged cor- 
relations correspond to weaker localization, ultimately al- 
lowing for critical delocalization at some particular value 
of 13. 

To conclude we have demonstrated a more physically 
transparent approach without the use of non-linear a 
models or renormalization group methods for analyz- 
ing the localization criteria for disordered tight-binding 
Hamiltonians with power-law hopping. To construct our 
perturbative approach we relied upon exactly summable 
series for the integer values of the hopping exponent a, 
allowing for a precise determination of the unperturbed 
dispersion relations. Then we utilized a simple scaling in- 
equality, which compares the fluctuations of the disorder 
to the energy level spacings of the unperturbed system. 
We have found that for a = 1, there are extended states 
and the perturbative condition holds for sufficiently weak 
disorder. 



As expected when ^q=4, which is more abrupt on the 
scale of the lattice, one recovers the scaling exponent of 
2/3 which is consistent with numerical simulations of dis- 
ordered ID chains^ with nearest neighbor hopping. It 
should be noted that the numerical prefactors depend 
on the specifics of the model and the proper choice of 
boundary conditions. 
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